Crossover from Thermal Activation to Quantum Interlayer Transport at the 
Superconducting Transition Temperature of Bi2Sr2CaCu208+5 
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We perform a detailed study of temperature, bias and doping dependance of interlayer transport 
in the layered high temperature superconductor Bi2Sr2CaCu208+i • We observe that the shape of 
interlayer characteristics in underdoped crystals exhibit a remarkable crossover at the supercon- 
ducting transition temperature: from thermal activation- type at T > Tc, to almost T— independent 
quantum tunneling- type, at T < Tc. Our data indicates that the interlayer transport mechanism 
may change with doping: from the conventional single quasiparticle tunneling in overdoped, to a 
progressively increasing Cooper pair contribution in underdoped crystals. 

PACS numbers: 74.72.Hs 74.45.-l-c 74.50.+r 74.25. Jb 



How does the high temperature superconductivity 
(HTSC) emerge with decreasing temperature? This 
highly debated question is crucial for understanding the 
mechanism of HTSC. The superconducting transition in 
HTSC appears to be unusual, with a seeming lack of 
changes in the quasiparticle (QP) density of states at 
Tc andpersistence of the normal state pseudogap above 
Tc [H, However, there is no consensus about the 
T— evolution of the energy gap. Results obtained by dif- 
ferent experimental techniques range from complete T— 
independence of the superconducting gap, A5G, in the 
whole T— range 3, to strong T— dependence at T = Tc 
d, H, H, 0, B 111- some experiments, coexistence of 
the two energy scales below Tc was reported 0, H, S ■ 
The existing controversy is one of the major obstacles for 
understanding HTSC and requires further T— dependent 
studies. 

Intrinsic tunneling spectroscopy utilizes the weak inter- 
layer (c— axis) coupling in layered HTSC. Unlike surface 
probe techniques, it is perfectly suited for T— dependent 
studies and provides information about bulk electronic 
properties of HTSC. In recent years different ways of im- 
proving this technique were employed 0, H, [Hf, [H, HI, 
[l4j . Yet, interpretation of intrinsic tunneling character- 
istics may not be free from misconceptions, associated 
both with the problem of self-heating [15] and the lack of 
complete understanding of the interlayer transport mech- 
anism. 

Here we study temperature, bias and doping dependen- 
cies of intrinsic tunneling in small Bi2Sr2CaCu208-i-i5 (Bi- 
2212) mesa structures. We observe an abrupt crossover 
in the shape of current- voltage characteristics (IVC's), 
from thermal activation (TA) like at T > Tc, to 
T— independent quantum tunneling like at T < Tc- Our 
data indicates that the interlayer transport mechanism 
in Bi-2212 changes in underdoped Bi-2212, together with 
development of the pseudogap. The observed simple TA 
behavior in the whole normal state region puts strong 



constrains on the nature of the c— axis pseudogap. In 
the supplementary information [l5| we discuss single QP 
characteristics, self-heating and clarify determination of 
parameters. 

We study two batches of crystals: the Y-doped Bi(Y)- 
2212 with the maximum Tc ~ 94. 5K and the pure Bi-2212 
with the maximum Tc ~ 86K. Details of mesa fabrication 
and characterization are described elsewhere 0, H, [13 ■ 

Figs. 1 a) and b) show the T— variation oi I — V and 
dl /dV{V) curves for a Bi-2212 mesa, containing N ^ 7 
intrinsic Josephson junctions, see Fig. 1 c). The follow- 
ing characteristic features are seen in dl/dV curves: The 
sharp "coherence" peak appears at T < Tc and is, in 
analogy with conventional low- Tc junctions, attributed 
to the sum-gap voltage Vg = 2NAsG/e 0, S, [O, [ll . A 
broader hump is seen in the whole T— range and is at- 
tributed to the c— axis pseudogap 0, d, [ll| . At elevated 
T the dl/dV curves cross in nearly one point [l,[l3]- The 
T— dependence of these features is plotted in Fig. 1 d) 

Figure 2 a) shows dI/dV{V) curves (in a semi- log 
scale) for the most underdoped Bi(Y)-2212 mesa with 
iV = 34 junctions. The thin and thick lines in Fig. 2 a) 
represent dl / dV {V) curves below and above Tc, respec- 
tively. It is seen that the shape of interlayer characteris- 
tics exhibits a remarkable qualitative change at Tc'. 

(i) Above Tc the slope of ln[dI/dV] (V) curves decreases 
with increasing T, while the voltage scale remains almost 
the same. The curves lean towards the horizontal line 
dl/dV — ctat, corresponding to the normal state conduc- 
tance, and cross almost in one point. 

(ii) Below Tc the ln[dI/dV]{V) curves have the same 
T — independent slope. With increasing T the curves re- 
main parallel and move towards the vertical axis V ^ 
(i.e., the voltage scale decreases). 

Experimental curves from Fig. 2 a) have a characteris- 
tic V-shape with voltage-independent slope. The overall 
shape of the curves at T > Tc closely resembles thermal 
activation characteristics [lB| , which can be described by 
a simple expression [l5l |: 
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FIG. 1: (Color online), a) /- V and b) dlldViy) characteris- 
tics for a Bi-2212 mesa at diflterent T . c) Periodic quasiparticle 
branches, d) Temperature dependence of the peak, hump and 
crossing voltages (per junction). 



where n is the concentration of mobile charge carriers, 
and Vta is the TA barrier. Indeed, the cosh term repro- 
duces the rounded V-shape of ln[c?//dF](y) curves with 
the slope that monotonously increases as l/T. 

In Fig. SIO [l^ we show that the data from Fig. 2 a) in 
the whole normal state, T > Tc, is very well described by 
Eq.([T]) with constant Uta — 2AmeV arid T— hnear n{T). 
It also reproduces the crossing point T6], which accord- 
ing to Eq.([T|) occurs at eV ~ 2Uta- Apparently, the 
TA barrier should be associated with the c— axis pseudo- 
gap. Thus, the appearance of the crossing point is simply 
the consequence of a fairly T— independent c— axis pseu- 
dogap y, 0, m ) while the crossing voltage represents 
the pseudogap energy. Indeed, a correlation between the 
hump and the crossing voltages is seen from Fig. 1 d). 

According to Eq.(IT]) the slope of the curves in Fig. 
2 a) should be proportional to the reciprocal temper- 
ature: d/dV{\n[dI/dV]) = (e/2fcsT) tanh(ey/2fcBT) ~ 
e/2kBT. Fig. 2 b) shows the effective temperature, 
Teff, obtained from the slopes at V/N =30mV 15]. It 
is seen that in the normal state T^ff ~ T, confirming 
the TA nature of interlayer transport at T > Tc- How- 
ever, at T < Tc an abrupt saturation of Teff (T) oc- 
curs, typical for quantum tunneling transport [3- 
This is the central observation of this work. Remark- 
ably, the crossover is clearly distinguishable at T — s- Tc, 
although no other spectroscopic features can be resolved 
in dI/dV{V) curves. 

So far we discussed dI/dV(V) at finite bias. Another 
straightforward way to investigate the TA behavior is 
to consider T— dependence of zero-bias resistance, Rq, 
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FIG. 2: (color online), a) A semi-log plot of dl /dV{V) curves 
for an underdoped Bi(Y)-2212 mesa at different T. Note that 
below Tc the curves maintain the same slope, while above 
Tc the slope changes progressively with T. A characteristic 
crossing point of dl /dV{V) curves at T > Tc is marked by the 
arrow, b) The effective temperature, obtained from the slopes 
of the dI/dV{V) curves at a finite bias V/N = 30mV [H. A 
remarkable crossover from the thermal activation, Tf.f f = T, 
to quantum, T^ff — const, behavior occurs at Tc. 



(the exp term in Eq.([T])). Experiments on small mesas 
provide a possibility to measure the small-bias QP re- 
sistance, R'^^ , in the superconducting state fill. figj . 
which is otherwise shunted by the supercurrent [20| . This 
is facilitated by the large specific capacitance (high qual- 
ity factor) of intrinsic junctions [l8[, which causes the 
pronounced hysteresis in the IVC's and allows junctions 
to remain in the resistive state even below the critical 
current, see Fig. 1 c). Therefore, bias yield an additional 
parameter for our studies, which may render crucial for 
correct interpretation of the data. 

Fig. 3 a) shows the dc~ resistance Rdc ^ V/I at differ- 
ent / for the same Bi-2212 meas as in Fig. 1. Measure- 
ments were done by first applying a large current, suffi- 
cient for switching to the last QP branch, see Fig. 1 c), 
and then ramping it back to the desired value. The Rdc 
at the smallest current / = 5fiA coincides with the con- 
ventional ac— resistance, Rac- It drops at T < Tc, since 
the finite quality factor of the junctions causes retrap- 
ping to the superconducting branch at a finite current 
[18|. From Fig. 3a) it is seen how non-linearity develops 
with decreasing both bias and T. At small / the curves 
approach an asymptotic, allowing a confident estimate 
of the zero-bias quasiparticle resistance, Rq^ , without a 
need for extrapolation. 
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FIG. 3: (Color online), a) T— dependence of the dc— resis- 
tance at different bias currents for the same Bi-2212 mesa. 
Panels b-d) show the asymptotic zero bias resistance: b) in 
the double-logarithmic scale; c) as a function of l/T; and d) 
normalized by T: a clear linear thermal activation behavior 
(dashed line) is observed in the whole normal state region. 
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FIG. 4: (Color online). T— normalized zero-bias quasiparticle 
resistance, extrapolated from IVC's at T < Tc and Raa 

measured with a small ac-current for Bi(Y)-2212 mesas with 
different doping. Dashed lines represent TA fits, Eq.((T}, at 
T > Tc- Inset shows relative values of the excess QP resis- 
tance with respect to the TA fit of the normal state. A large 
excess resistance appears below Tc in the overdoped mesa. 
The excess resistance rapidly decreases with decreasing dop- 
ing and becomes negative in the underdoped mesa. 



The R^"^ grows with decreasing T at T < and 
is usually described either in terms of power-law depen- 
dence [20], inherent for single QP tunneling in the pres- 
ence of a d— wave gap (see Fig. S7 [3]), or in terms of 
TA The double- logarithmic plot, Fig. 3 b), demon- 
strates that Ro{T) is not described by the power law in 
any extended T— range for our mesas. Neither is it per- 
fectly described by the Arrhenius law explUTA/ksT), 
as demonstrated in Fig. 3 c). On the other hand. Fig. 
3 d) demonstrates that the ratio Rq/T follows very ac- 
curately the Arrhenius law (dashed line) in the whole 
normal region, consistent with the TA expression Eq.([T]) 
and with the finite bias behavior. Fig. 2 b). However, at 
T < Tc, Rq^ /T deviates downwards from the Arrhenius 
law. This is consistent with saturation of the effective 
Teff{T < Tc), as shown in Fig. 2 b). Therefore, it is a 
consequence of the same crossover. 

In Fig. 4 we analyze doping dependence of Ro{T)/T. 
In the normal state, Ro{T) varies in the TA manner for 
all doping levels. In the superconducting state, Rq^ also 
continues to grow with decreasing T, but the rate of the 
growth with respect to the TA behavior, Rta, (dashed 
lines) strongly depends on doping: 

For the overdoped (OD) mesa, Rq^{T < Tc) grows 
much faster than in the normal state, i.e., there is a sharp 
onset of the excess resistance at T = Tc- However, the 
excess resistance at T < T^, decreases rapidly with un- 
derdoping. It becomes small in the near optimally doped 



(OP) mesa and turns negative already in moderately un- 
derdoped (UD) mesas. 

Inset in Fig. 4 shows the magnitude of the excess QP 
resistance with respect to the TA fit, Rta, in the nor- 
mal state (dashed lines). The quality of the fit could 
be judged from the flatness of the normal state region. 
A progressive decrease of the excess QP resistance at 
T < Tc with decreasing doping is obvious. 

We start discussion of the observed phenomena by rec- 
ollecting the expected T— variation for conventional sin- 
gle QP tunneling characteristics. As discussed in sec. I 
of the supplementary information, provided the physi- 
cal requirement that Asg and the quasiparticle lifetime 
increase with decreasing T, all single QP tunneling char- 
acteristics (including d— wave) exhibit the following uni- 
versal T— dependence: 

A) Opening of the superconducting gap at T < Tc 
leads to appearance of the large excess QP resistance due 
to rapid freezing out of QP's upon their condensation into 
Cooper pairs (see inset in Fig. S7 p^). 

B) Rq^ continues to grow with decreasing T both due 
to freezing-out of thermal QP's and due to growth of the 
QP lifetime, which decreases the number of available sub- 
gap QP states. All this leads to a progressive growth of 
the slope of dI/dV{V) characteristics with decreasing T. 

The sudden appearance of excess resistance at T < Tc 
in OD mesas indicates that the superconducting transi- 
tion here is conventional, in a sense that it is accompanied 



4 



by opening of the superconducting gap, and is also con- 
sistent with single QP tunneling mechanism of interlayer 
transport. 

However, the superconducting transition in under- 
doped mesas is clearly abnormal. We emphasize again 
that appearance of the negative excess resistance and the 
crossover to T— independent slope reflect the same phe- 
nomenon at zero and finite bias, respectively. Therefore, 
any explanation of the observed unusual T— dependence 
must account for both zero and finite bias behavior. 

At the first glance, the lack of excess resistance at 
T < Tc in UD mesas could be consistent with the pre- 
cursor superconductivity scenario of the pseudogap [sj, 
according to which the gap does not open and QP's do 
not start to pair at Tc, but at much higher T*. Then 
one could argue that the lack of excess resistance reflects 
simply the lack of dramatic changes in the QP spectrum 
upon establishing of the global phase coherence at Tc- 
However, the same argument would make it difficult to 
explain the abrupt crossover to T— independent slope. 
Furthermore, as we noted above, the slope of dI/dV{V) 
should continue to increase with decreasing T for single 
QP tunneling. Therefore, the observed crossover to con- 
stant slope can hardly be explained in terms of single QP 
tunneling [15j . 

This conclusion brings us to one possible interpreta- 
tion of the observed phenomena. If the c— axis transport 
in UD mesas is not solely due to single QP's, it must also 
involve pairs. The corresponding multi-particle tunnel- 
ing [22, the multiple Andreev reflection [H, [2^ 
processes are well studied for conventional low— Tc junc- 
tions. Both processes are similar [25!] and are almost T— 
independent, because they do not rely on the presence of 
thermally excited QP's, but directly involve Cooper pairs 
from the Fermi level. The multi-particle current occurs 
via elastic conversion (dissociation or recombination) of 
Cooper pairs into QP's; while the multiple Andreev re- 
flection process is due to transmutation of a quasielectron 
into a quasihole with creation of a Cooper pair and can 
be universally described in terms of inelastic tunneling in 
the presence of time- dependent phase difference due to 
the ac- Josephson effect [25:] . 

Importantly, the multi-particle processes result in al- 
most T— independent and F— exponential subgap current 
[H, US] , which essentially coincides with the description 
of the observed crossover, see Fig. 2. Furthermore, the 
abruptness of the crossover at Tc, points towards the 
coherent multiple Andreev reflection mechanism of the 
interlayer current, because it requires the time-periodic 
ac- Josephson effect [25| and, consequently, abruptly dis- 
appears simultaneously with phase coherence at Tc- 

We may further speculate why the multiparticle pro- 
cesses become progressively more important with under- 
doping. The multiparticle current decrease much faster 
than the single QP current with decreasing the interface 
transparency [IJ, [2^. Therefore, the single QP domi- 
nates over multiparticle current, unless there are micro- 
shots in the tunnel barrier [23[. The required micro- 



shots in underdoped Bi-2212 could arise from doping- 
dependent nano-scale inhomogeneity, observed at the 
surface of Bi-2212 (26l[27|. which in this case should per- 
sists also in the bulk of Bi-2212 crystals. 

Such interpretation implies that not only the elec- 
tronic structure, but also the c— axis transport mechanism 
changes with underdoping: from coherent and directional 
[m single QP tunneling in overdoped, to progressively 
increasing pair contribution in underdoped Bi-2212. The 
latter is apparent only in the phase coherent state at 
T < Tc, is almost T— independent, and is consistent with 
multiple Andreev reflection mechanism of the interlayer 
transport. 

Finally we note that the reported amazingly trivial 
thermal activation behavior in the whole normal state re- 
gion, puts strong constrains on the nature of the c— axis 
pseudogap. Even though, the similar TA-behavior can 
be ascribed to several processes (e.g., inelastic tunneling 
via the impurity ^27i . i28] , or elastic tunneling via a reso- 
nant [2^ state in the tunnel barrier, or Coulomb blocking 
of tunneling it does not involve any gap, nor angu- 
lar dependence in the QP spectrum, but assumes instead 
that there is some constant blocking barrier for interlayer 
hoping [1^. This may indicate that the "large" c— axis 
pseudogap is not a pairing gap, which would naturally ex- 
plain it's indifference to such classical depairing factors 
as temperature and magnetic field [7,]. 

We are grateful to E.Silva for stimulating discussions 
and to A.Yurgens, T.Benseman and G.Balakrishnan for 
providing Bi(Y)-2212 and Bi-2212 single crystals. Finan- 
cial support from the K.&A. Wallenberg foundation and 
the Swedish Research Council is gratefully acknowledged. 

I. SUPPLEMENTARY INFORMATION 

A. Nomenclature 

ARPES: angular resolved photoemission spectroscopy, 

Bi-2212: BizSraCaCuaOg+s, 

Bi(Y)-2212: Y-doped Bi-2212, 

DoS: density of states, 

A: energy gap, 

AT: self-heating, 

F: depairing factor, reciprocal quasiparticle lifetime; 
HTSC: high temperature superconductor, 
J: current, 

IJJ: intrinsic Josephson junction. 
IVC: Current- Voltage characteristic, 
N: number of IJJ's in the mesa, 
n: concentration of mobile charge carriers, 
CD, OP, UD: overdoped, optimally doped, under- 
doped, respectively; 
QP: quasiparticle, 

Rq: zero bias resistance (in the normal state), 
Rq : zero bias quasiparticle resistance in the super- 
conducting state, is different from Rq because of the pro- 
nounced hysteresis in the IVC's at T < Tc; 
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i?„: tunnel (normal) resistance, 

Rth- thermal resistance of a mesa, 

SIS: Superconductor-Insulator-Superconductor, 

SIN: Superconductor-Insulator-Normal metal, 

STM: scanning tunneling microscope, 

TA: thermal activation, 

Tc- superconducting critical temperature, 

Uta- thermal activation barrier, 

V: voltage, 

Vg-. the sum-gap voltage Vg — 2NA/e 

The capital M/S in front of the figure or equation num- 
ber indicates that the reference is made to the Manuscript 
or the Supplementary information, respectively. 



B. Influence of the quasiparticle life time on dl/dV 
characteristics 

We have argued in the manuscript that the observed 
crossover to T— independent slope and the correlated ap- 
pearance of the negative excess resistance are inconsistent 
with the single QP tunneling mechanism of the c— axis 
transport. However, the interlayer tunnelin g in Bi-2212 
may depend on a number of parameters [30. l3ll. [3^ : 

(i) The QP life time, l/F, which, according to ARPES 
data, has a substantial T— dependence close to Tc, where 
it increases roughly linearly with T; [s^ 

(ii) Momentum conservation upon tunneling (coherent 
vs. incoherent tunneling); 

(iii) The directionality of c— axis tunneling, i.e., the 
angular dependence of the tunneling matrix element 
i^i,(/32), caused by a non-spherical Fermi surface (sil . 

(iv) Temperature and (v) angular dependence of the 
gap in the DoS and the points (i-iii) above. Deviations 
from pure d— wave symmetry were reported in recent 
ARPES experiments %. 

It may not be obvious that experimental data can not 
be described by the single QP tunneling with a fortunate 
combination of all those parameters. Therefore, below 
we want to dwell upon this statement. 

We start with considering the consequences of points 
(ii, iii) and (v), i.e., coherence and directionality of tun- 
neling and the symmetry of the order parameter. The 
single QP tunneling current is given by: 

/'27r /'27r /' + cxd 

I{V)= / / / d^id^2dE (2) 

Jo Jo J-oo 

t{^i, f2)p{E, fi)p{E + eV, V2)I{E) [1 - J{E + eV)] , 

where E is the energy of the QP with respect to the 
Fermi surface, 1^91^2 are the angles in the momentum space 
of the initial and the final state of the QP, and p{E, ipi) 
p{E + eV, LP2) are the corresponding QP DoS: 

p{E,v)^^[{E-iT)/^{E~^W~X{^l (3) 



In Fig. S5 a) and b) we show numerically simulated 
SIS characteristics for different single QP tunneling sce- 
narios. We assumed, t{ifi^if2) —const for incoherent- 
nondirectional, t{ipi,ip2) oc 5{Lpi — (^2) for coherent- 
nondirectional and t((pi, (^2) oc [cos(fca;) — cos(fcy)]^(5((^i — 
ip2) for coherent-directional tunneling [3l|. The simu- 
lations were made for A((p = 0) = 35meV, at low T 
and small depairing F <C A(0). Detailed discussion of 
dl/dV characteristics can be found in Ref. [30|. Irre- 
spective of the scenario, the zero-bias resistance diverges 
(the conductance tends to zero) at T — > and F = 0. 
The divergence is removed by the finite F, see the inset 
in Fig. S6 b). Therefore, dI/dV{V = 0,T ^ 0) pro- 
vides the information about the value of F. We also note 
that dI/dV{V) characteristics exhibit a sharp peak at 
the sum-gap voltage Vg = 2A/e ( except for the curve 
A). 

It is expected that interlayer QP tunneling in Bi-2212 
single crystals should be predominantly coherent (pro- 
vided that the single crystal is pure enough so that there 
is no momentum scattering upon tunneling) and strongly 
directional, with dominating antinodal tunneling ^1, 331. 
Indeed, the curve C in Fig. S5 resembles most closely 
the experimental characteristics (cf. with the curves at 
T = 13 AK from Fig. Mia and b), see also a discussion 
inRefs. 

In what follows, we will analyze the incoherent- 
nondirectional single QP tunneling characteristics, rep- 
resenting an "extreme" d— wave case with the weakest 
sum-gap peak. The conclusions will, however, be univer- 
sal for any single QP tunneling scenario. 

Now we consider the effect of T (point iv) and F (point 
i) on dl/dV characteristics. It is clear that both param- 
eters will smoothen dI/dV{V) characteristics and fill-in 
the dip in conductance at y = 0. However, they do this 
in a slightly different manner. This is illustrated in Figs. 
S6 and S7. 

The increase of F leads to smearing of the QP gap sin- 
gularity and simultaneous filling-in of the sub-gap states 
in DoS, Eq.(S2). This naturally leads to smearing of 
the sum-gap peak and to filling-in of the zero-bias dip in 
dl/dV, as shown in the inset in Fig. S6 b). 

The increase of T leads to the decrease of A. Both fac- 
tors result in the increased number of excited QP's above 
the gap, which leads to rapid filling-in of the zero bias 
dip in conductance. However, the increase of T (unlike 
F) keeps the gap singularity in DoS unchanged. From 
Fig. S7 it is seen that this leads not only to filling-in 
of the dip, but also to development of the maximum at 
V = 0, representing the so called zero-bias logarithmic 
singularity (35j . 

The origin of the singularity is very simple: at elevated 
T there is a substantial amount of thermally excited QP's 
just above the gap. At V = the partly filled gap sin- 
gularities in the two electrodes are co-aligned, causing a 
large current flow from one electrode to another, which is 
exactly compensated by the counterflow from the second 
electrode. However, exact cancelation is lifted at an arbi- 
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FIG. 5: (Color online). Simulated single QP characteris- 
tics for different tunneling scenarios. A: incoherent, non- 
directional, d— wave; B: coherent, non-directional, d— wave; 
C: coherent, directional, d— wave; D: s— wave. Note that the 
sharp sum-gap peak and the half-gap singularity in dl/dV are 
inherent to most scenarios. 



FIG. 6: (color online). Effect of the depairing factor on single 
QP characteristics. Shown are simulated characteristics at 
constant T and A for varying F. Increase of F tends to reduce 
the slope of dl/dV characteristics. Inset in panel b) indicates 
that filling-in of the zero bias dip with F occurs in a power 
law manner. 



trary small voltage across the junction, leading to a sharp 
maximum in dl/dV. From Fig. S7 it is seen that the 
zero-bias maximum in d— wave junctions is pronounced 
even in the extreme case of weakest sum-gap peak. We 
have checked that the zero-bias logarithmic singularity is 
inherent for single QP characteristics, irrespective of the 
tunneling scenario. 

Interestingly, the experimental characteristics do not 
exhibit the zero-bias logarithmic singularity, see Fig. 
M2a). Within the single QP tunneling scenario this is 
only possible if the depairing factor increases with T at 
such a rate that it smears out the gap singularity in DoS 
and thus suppresses development of the singularity. We 
estimate that the corresponding F at T = Tc should be 
in the range 2-5 meV, somewhat smaller than deduced 
from ARPES [i,[3|. 

Now we can substantiate our statement. In Fig. S8 
we show the attempt of maintaining the same slope of 
lndI/dV(V) curves for progressively decreasing A. Since 
increase of both F and T leads to filling-in of the zero bias 
dip, i.e., to decrease of the slope, the slope can be main- 
tained only if F and T move in the opposite directions. 



It is impossible to maintain the slope under the physical 
requirements that F increases and A decrease simulta- 
neously with increasing T because all of the three pa- 
rameters work in the same direction - decrease the slope. 
The inset in Fig. S7 shows the calculated zero bias re- 
sistance vs. 1/T for incoherent, non-directional, d— wave 
single QP tunneling and different F. It is seen that a 
prominent excess resistance appears even in this extreme 
d— wave case and even for very large F. The dashed line 
shows that the excess resistance grows approximately as 
T^^ in this case. The abrupt appearance of the excess 
QP resistance is a mere reflection of the abrupt opening 
of the superconducting gap at Tc- Therefore, the ob- 
served T— independent slope and the absence of excess 
resistance, observed for UD mesas, are inconsistent with 
the single QP tunneling and points towards the doping- 
dependent change in the interlayer transport mechanism. 
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FIG. 7: (color online). Effect of T and A(r) variation on 
single QP characteristics. Calculations were made for inco- 
herent, non-directional, d— wave tunneling, the experimental 
A(r) dependence and linear T{T). Note that increase of T 
and decrease of A leads not only to smearing of dl /dV charac- 
teristics, but also the appearance of the logarithmic zero-bias 
singularity. Inset shows the zero bias resistance vs. l/T for 
different values of r(rc). It is seen that the excess QP resis- 
tance rapidly develops at T < Tc even in the d— wave case. 




Incoherent, 
non-directional 
d-wave tunneling. 
r=1, A=35 (meV), T= 59K; 
r=2, A=35 (meV), r=55.2K; 
r=5, A=35 (meV), T= 36. 2K; 
=5, A=17.5 (meV), T= 32.4K; 
=7, A=10(meV), r=0.1K. 
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FIG. 8: (Color online). An attempt to maintain the slope of 
single QP characteristics at different T. It is seen that the 
slope can be maintained only under the un-physical require- 
ment that r increases and A decreases with decreasing T. 



C. Analysis of self- heating and non-equilibrium 
phenomena 

Self-heating can distort IVC's of Josephson junctions 
at large bias. The temperature rise due to self-heating is 
given by a simple expression 

^T = PRth{T), 



where P = IV is the dissipated power and Rth is the 
effective thermal resistance of the junctions, which is 
T— dependent and, therefore, bias dependent [11]. The 
influence of self-heating on IVC's of our mesa structures 
was thoroughly studied in Rcfs.fsl, [H, [13, [11] (see also 
a discussion of scaling of IVC's in the inset of Fig. 1 and 
Fig. 3 from Ref.[7^). 

Despite relative simplicity of the phenomenon, discus- 
sion of self-heating in intrinsic tunneling spectroscopy has 
caused a considerable confusion, a large part of which has 
been caused by a series of publications by V.Zavaritsky 
[3^ , in which he " explained" the non-linearity of intrin- 
sic tunneling characteristics by assuming that there is no 
intrinsic tunneling. The irrelevance of this model to our 
subject was discussed in Ref . (sst . 

A certain confusion might be also caused by a large 
spread in Rth, reported by different groups [3, [IE [iOj 
[ill . |43 |. For the sake of clarity it should be empha- 
sized that those measurements were made on samples 
of different geometries. It is clear that Rth depends 
strongly on the geometry [H, [13|, e.g., Rth can be 
much larger in suspended junctions with poor thermal 
link to the substrate [i^ than in the case when both 
top and bottom surfaces of the junctions are well ther- 
mally anchored to the heat bath [l3|. For mesa struc- 
tures similar to those used in this study (a few /xm in- 
plane size, containing iV ~ 10 IJJ), there is a consensus 
that Rth{4:.2K) - 30 - 7QK/mW (depending on bias) 
[lE[4l| and Rth{90K) ~ 5 - lOK/mW [11]. Larger val- 
ues Rth > lOOK/mW claimed by some authors [lo] are 
unrealistic for our mesas because they can withstand dis- 
sipated powers in excess of lOmW without being melted. 

Yet, we note that talking about a typical value of Rth is 
equally senseless as talking about a typical value of a con- 
tact (Maxwell) electrical resistance: both depend on the 
geometry. Therefore, reduction of mesa sizes provides 
a simple way for reduction of self- heating [13]. Conse- 
quently, variation of dl /dV characteristics with the junc- 
tion size and geometry provides an unambiguous way of 
discriminating artifacts of self-heating from the spectro- 
scopic features (361 ]. 

How self- heating can distort the IVC's of Josephson 
jimctions is obvious: since self- heating rises the effective 
T it may affect the IVC only via T— dependent parame- 
ters. There are three such parameters: t he q uasiparticle 
resistance, the superconducting switching 18, 'is^ current 
and the superconducting gap. They will affect the IVC's 
of Bi-2212 mesas, containing several stacked Josephson 
junctions, in the following manner: 

The consecutive increase of T upon sequential switch- 
ing of IJJ's from the superconducting to the resis- 
tive state will distort the periodicity of quasiparticle 
branches. Each consecutive QP branch will have a 
smaller QP resistance (smaller V at given I) and smaller 
switching current, see the discussion in Refs. [3,[lHl- This 
type of distortion of the QP branches becomes clearly 
visible (at base T = ^.2K) when AT > 2QK 

The T— dependence of A may lead to appearance of 
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FIG. 9: (Color online). Simulated distortion of SIS tunneling 
characteristics by strong self-heating, a) Undistorted IVC's 
at different T. Simulations were made for typical parameters 
of our mesas, T— dependent thermal conductivity of Bi-2212, 
and for coherent, directional, d— wave tunneling, b) distorted 
IVC's at the same base temperatures; c) The mesa tempera- 
ture as a function of bias, d) Temperature dependence of the 
genuine superconducting gap (solid line) and the "measured" 
gap obtained from distorted IVC's (dashed curve). Note that 
even strong self-heating (T reaches Tc/2 at Vg at 4.2K) does 
not cause considerable distortion of the measured gap. Data 
from Ref.fsl 



back-bending of the IVC at the sum-gap knee. Fig. S9 
reproduces the results of numerical simulation of such 
the distortion, made specifically for the case of Bi-2212 
mesa with the corresponding T— dependent parameters 
(see Ref. j34j for details). Fig. S9 a) shows a set of simu- 
lated undistorted IVC's at different T for coherent, direc- 
tional, d— wave tunneling with a A(T), shown by the solid 
line in Fig. S9 d). Panels b) and c) show the distorted 
IVC's and the actual junction temperature, respectively. 
The dashed line in panel d) represents the "measured" 
gap obtained from the peak in distorted dl /dV character- 
istics. Remarkably, the deviation from the true A(T) is 
marginal, despite large self-heating, AT ~ Tc/2 at 4.2X1 

The robustness of the measured gap with respect to 
self-heating is caused by a flat T— dependence of the su- 
perconducting gap at r < Tc/2. If self-heating reaches 
Tc at the sum-gap knee, an acute back-bending appears 
in the IVC's, however even this does not cause princi- 
ple changes in the "measured" A(T). In experiments on 
large mesas (i^l or suspended structures [i^, in which 
acute self-heating was reported, no clear Ohmic tunnel- 
ing resistance could be observed at high bias, in stark 
contrast to our mesas, see Fig. Mia, and Refs.0, H, 

Simulations as in Fig. S9 clearly show that irrespective 



of self-heating the "measured" gap vanishes at T = Tc, 
i.e., simultaneously with the true gap. Thus trivial self- 
heating simply can not "hide" the qualitative A(T) de- 
pendence. For example, there is no way in which one can 
get the vanishing "measured" (self- heating affected) gap 
if the true gap is T— independent. Therefore, the discrep- 
ancy in the measured strong T— dependence of the gap 
at Tc in intrinsic tunneling (and recent ARPES [9]), or 
complete T— independence in STM measurements [3| can 
not be attributed to self-heating. 

To quantitatively analyze the significance of self- 
heating in our mesas, we provide values of the dissipation 
power P = IV for the studied mesas: For the near-OP 
Bi(Y)-2212 mesa at T = A.QK dissipation at the sum- 
gap peak is 0.82mW. From the previous analysis of 
the size dependence of Vg for similar mesas [36[ it was 
observed that Vg becomes size-independent, and there- 
fore not affected by self-heating, for small mesas with 
P{Vg) < ImW. All the mesas studied in this manuscript 
fall into this category. All of them exhibit perfect pe- 
riodicity of quasiparticle branches fi\ and none of them 
exhibit back-bending at any T, which according to Fig. 
S9, implies that self-heating at Vg is less than Tc/2 at 
T = A.2K. This is consistent with previous in-situ mea- 
surements dil, Rth{T = A.2K, F = I/g) ~ iQ-mK/mW. 

The rest of the data presented in the manuscript is not 
affected by self-heating for the following reasons: 

The effective Tg// in Fig. M2b) was obtained at 
V/N — iQmeV at which the dissipation power was only 
P ~ i.9,^lW at T = 5.5if (for comparison, P{Vg) = 
Q.21mW for the same curve). Therefore, self-heating 
here is negligible (sub-Kelvin). Besides, as argued in the 
manuscript, the slope of hidl /dV{V) in Fig. M2 is ap- 
parently V and P independent, therefore the same data 
could be obtained at much smaller or larger P. The data 
in Figs. M3 b-d) and M4 is free from self-heating because 
it was obtained at zero bias. 

Finally, it is important to emphasize that the concept 
of heat diffusion is inapplicable for small Bi-2212 mesas 
containing only few atomic layers. The phonon trans- 
port in this case is ballistic [la, EH] and the energy fiow 
from the mesa is determined not by collisions between 
the tunneled non-equilibrium QP with thermal phonons, 
but by spontaneous emission of a phonon upon relax- 
ation of the non-equilibrium QP [13]. This process is not 
hindered at T 0. Therefore, the effective Rtt (and self- 
hcating) can be much smaller because it is not limited by 
poor thermal conductivity at T = 0, but is determined 
by the fast, almost T— independent, non-equilibrium QP 
relaxation time. The concept of self-heating becomes ad- 
equate only in the bulk of the Bi-2212 crystal, where the 
dissipation power density and the temperature rise are 
much smaller due to the much larger area of the crystal. 
For more details see the discussion in Ref. [s^ . The non- 
equilibrium energy transfer channel is specific for atomic 
scale intrinsic Josephson junctions made of perfect single 
crystals. It can explain a remarkably low self-heating at 
very high bias p^ . 
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FIG. 10: (Color online), a) Simulated dl/dV characteristics 
according to Eq.(Ml) for the same T as in Fig. M2 a). It is 
seen that at T > Tc they provide a good fit to experimental 
data, including the crossing point, b) The slope of experi- 
mental characteristics from Fig. M2 a) at V/N = 30mV (red 
squares, left axis) and the logarithm of zero bias resistance 
(green circles, right axis) vs. 1/T. It is seen that both fol- 
low the TA behavior at T > Tc and simultaneously deviate 
downwards in the superconducting state, indicating that ap- 
pearance of the negative excess resistance and the crossover to 
quantum c— axis transport are correlated. The dashed-dotted 
line represents the slope of resonant tunneling characteristics 
(shown for comparison). 



D. Thermal activation analysis 

The TA current is given by a simple expression: 

Uta 



It A oc n(T) exp 



cosh 



eV 



2kBT 



(4) 



from which follows Eq.(Ml). 

Fig. SIO a) shows the dI/dV{V) characteristics in a 
semi-log scale calculated from Eq.(Ml) for the same T 
and scale as in Fig. M2 a), in order to facilitate the di- 
rect comparison. Simulations were made for constant ctjv 
and Uta = 2AmeV , which were adjusted to obtain sim- 
ilar Ro{T) values as in the experimental data from Fig. 
M2 a) at T > Tc; and for n{T) oc T. Considerable grows 
of mobile carrier concentration n{T) with T in UD crys- 
tals was reported in Hall effect measurements [i^. The 
assumption of linear n{T) is supported by observation 



that for strongly UD mesas it is Rq (circles in Fig. SIO 
b), rather than Rq/T, that is more accurately described 
by the Arrhenius law (dashed line). 

From comparison of Figs. M2a) and SlOa) it is seen 
that Eq.(Ml) does provide a good fit to experimental 
data at T > T^, including the crossing point. Simulta- 
neously a dramatic discrepancy is seen below Tc- This 
indicates that the slopes of dI/dV{V) curves in the su- 
perconducting state are not determined by the real T, 
but by some constant effective temperature, as shown in 
Fig. M2 b). 

Fig. SIO b) shows the slope of experimental curves 
from Fig. M2a) at V/N = 30mU (red squares) and 
the logarithm of zero bias resistance (green circles) as 
a function of 1/T. Clear linear dependence of both at 
T > Tc indicates that the IVC's in the normal state 
have thermal activation nature both at zero and finite 
bias, consistent with Eq.(Ml). According to Eq.(Ml) 
the effective temperature can be explicitly obtained from 
the slope d/dV[\n{dI /dV)] = (e/2fcBT) tanh(ey/2fcBT) 
~ e/2kBT. Here there are two slight obstacles: First, 
the tanh term deviated from unity at high T leading 
to eventual saturation of the slope at 1/T 0. The 
exact T— depensence according to Eq.(Ml) is shown by 
the solid (magenta) line in Fig. SIO b). We emphasize 
that there are no fitting parameters in this curve. It is 
seen that Eq.(Ml) reproduces quite well T— dependence 
of the slopes, except for a small offset. This additional 
offset is caused by the fact that the slope of experimen- 
tal curves becomes negative at large T Q. The simple 
TA expression, Eq.(Ml), does not reproduce this nega- 
tive slope. However, more advanced TA simulations do 
reproduce the negative slope and, in fact, all the details 
of experimental dl/dV characteristics at T > Tc [iGj . 

It should be noted that TA-like behavior is quite uni- 
versal for many process. Except for pure thermal ac- 
tivation over the finite barrier without the gap in the 
electronic DoS [T6| . it may appear in pure tunneling 
characteristics in the presence of the gap in electronic 
DoS, due to TA-like behavior of the Fermi factor (even 
though this would require a very specific correlation be- 
tween T— dependent factors mentioned in sec. I); as a re- 
sult of inelastic tunneling via the impurity [2^ , or elastic 
tunneling via a resonant state [2^ in the tunnel barrier. 
Abrikosov has shown that the later scenario can quanti- 
tatively reproduce the interlayer characteristics of HTSC 
in the normal state [2^ . The blue dash-dotted line in Fig. 
SIO b) represents the slope d/dV[ln{dI /dV)] in the case 
of resonant tunneling with the appropriate energy of the 
resonant state. It follows the simple linear TA behavior 
at ksT lower than the resonant energy and saturates at 
higher T. Comparison with the resonant tunneling cal- 
culation shows that the saturation of the slope at high 
T results in appearance of a negative offset in the linear 
d/dV[ln{dI/dV)] vs 1/T dependence at 1/T ^ 0. Be- 
low the saturation temperature, the actual temperature 
can be easily extracted from the d/ dV[ln{dI / dV)] slopes, 
compensated by the offset at 1/T — > 0, as indicated in 
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Fig. SIO b) : 



2kB [d/dV[ln{dI/dV)] - Offset{T-^ 0)] ' 

This expression was used for obtaining the effective T 
shown in Fig. M2b). 

From Fig. SIO b) it is seen that at T < Tc both lni?o 



and the slope of experimental curves simultaneously de- 
viate downwards with respect to the linear TA behav- 
ior. Therefore, both are consequences of the same phe- 
nomenon (which we attributed to doping-induced change 
in the c— axis transport) only at zero and finite bias, re- 
spectively. 
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